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Employing virtual crystal approximation and super-cell methods for doping, we have performed a 
comparative study of the electronic structures of various doped BaFe 2 As 2 materials by first principles 
simulations. Both of these methods give rise to a similar density of states and band structures in 
case of hole doping (K doping in Ba site) and iso-electronic P doping in As site. But in case of 
electron doped systems with higher doping concentration, electronic structures, calculated using 
virtual crystal approximation approach deviates from that of the super-cell method. On the other 
hand in case of iso-electronic Ru doping implemented by virtual crystal approximation, an extra shift 
of the chemical potential in electronic structure in comparison to super-cell method is observed and 
that shift can be used to predict the correct electronic structure within virtual crystal approximation 
as reflected in our calculated Fermi surfaces. But for higher Ru doping concentration, simple 
shifting of chemical potential does not work as the electronic structure calculated by virtual crystal 
approximation approach is entirely different from that of the calculated by super-cell formalism. 


I. INTRODUCTION 

Advent of high temperature superconductivi^ in Fe- 
based compounds in the vicinity of magnetism [iMSl has 
led to immense intrigue among the researchers to resolve 
the unsolved mysteries in the field of high temperature 
unconventional superconductivity. Apart from some old 
issues like superconducting pairing symmetry, proxim¬ 
ity of superconductivity and magnetism, unusual Fermi 
surface (FS) topology, a number of new phenomena like 
existence of orbital density wave, origin of nematic phase 
etc., are still under substantial debate [Ml. In liter¬ 
ature there are large number of experimental as well as 
theoretical studies regarding the electronic structures and 
properties of Fe-based superconductors [illiMl to es¬ 
tablish the mechanism of unconventional superconductiv¬ 
ity with high Tc. It should be mentioned that the highest 
transition temperature (Tc) till now for Fe-based super¬ 
conductors is 109K [ 2 ^. Correlation of superconduct¬ 
ing transition temperature and other physical properties 
with structural parameters is well established ( 2 ll - l^ . 
Fe-As(pnictides)or Fe-Se(chalcogenides) planes are back 
bone structures of these superconductors where As atoms 
are placed above and below the Fe plane. Height of this 
As atom from Fe plane is defined as ’anion height’ closely 
related to zas i.e., the fractional co-ordinate of As atom. 
It is also evident from literature that, as far as veracity of 
calculation of electronic structure is concerned, the role of 
structural parameters specially zas is instrumental (25I - 
. Fe-based superconductors are broadly classified into 
six families among them 122 family is the most studied 
family, both experimentally as well as theoretically. The 
main reason for that is the availability of high quality sin¬ 
gle crystal for 122 systems. It should also be noted that 
the parent compound BaFe 2 As 2 is not superconducting 
in ambient pressure but it shows magnetic order in the 


form of spin density wave (SDW). Superconductivity ap¬ 
pears in this system when SDW order is suppressed and 
that can be done either by applying hydrostatic pressure 
[ 2 ^ 4 ^ or chemical pressure ie., by doping. Supercon¬ 
ductivity can arise into the system by doping on any of 
the three sites including Fe sites Another in¬ 

teresting case is isoelectronic Ru substitution in place of 
Fe which also forge the system to become a supercon¬ 
ductor [13, ■ Electronic structure of all these doped 

122 systems has been studied using first principles calcu¬ 
lation based on density functional theory (DFT) with a 
great success [I 3 , @ . But there are some challenges 

regarding the application of DFT in Fe-based supercon¬ 
ductors whose electronic structure is very much sensitive 
to temperature, pressure and doping. We extend this is¬ 
sue in the theory and computational method section in 
a greater detail. 

As mentioned earlier in BaFe 2 As 2 system, one can 
dope on any of the sites to make it superconducting. 
Phase diagrams of Fe-based superconductors play a cru¬ 
cial role in the understanding of various physical prop¬ 
erties, which are very sensitive to doping and tempera¬ 
ture. A generic phase diagram of 122 family of Fe-based 
superconductors consists of various phases like supercon¬ 
ductivity, SDW order, C 2 and C 4 structural transition, 
nematic phase etc. (39l-l4^. All these phases are inti¬ 
mately dependant on doping concentration. There exists 
a large diversity in the phase diagrams of various doped 
BaFe 2 As 2 systems. For electron doped system supercon¬ 
ducting and other exotic phases exist in the system upto 
15-20% doping concentration which is very less compared 
to other 122 systems [43 - [43| . Theoretically this dop¬ 
ing effect can be accomplished by several ways but two 
of them are widely used: super-cell (SC) method and 
virtual crystal approximation (VGA). We later discuss 
these two methods concisely. Application of these two 
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methods to implement doping for calculating electronic 
structures (density of states, band structures, Fermi sur¬ 
faces etc.) of various doped 122 systems are there in the 
literature. But a comparative study of electronic struc¬ 
tures calculated using both these two approaches is rare. 
In this work we proclaim the differences in electronic 
structures of the various doped 122 systems calculated 
by implementing doping in VGA and SC methods. Five 
doped BaFe2As2 systems viz. Bai_2:K2;Fe2As2 (hole dop¬ 
ing), BaFe2-xCo3,As2 (electron doping), BaFe2-a;Nia;As2 
(electron doping), BaFe2(Asi_xPx)2 (isovalent doping) 
and BaFe2-a;Rua;As2 (isovalent doping) has been exam¬ 
ined through detail study of density of states and band 
structures using VGA and SC approaches. In particular, 
we observe that both these methods give rise to very dif¬ 
ferent electronic structures when one dope on Fe site. We 
also found that for higher doping concentration electronic 
structures calculated using VGA method deviate from 
that of the calculated utilizing SC formalism for certain 
cases. Our work demonstrate a detailed study of elec¬ 
tronic structures of all these above mentioned doped 122 
systems where doping has been implemented by two dif¬ 
ferent theoretical formalisms. In order to assert a correct 
electronic structures of these doped systems one must as¬ 
sure the validity of VGA formalism and one should also 
be aware of the limitation of these two methods. 
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FIG. 1: (Colour online) Calculated orbital projected partial 
density of states (PDOS) of undoped BaFe 2 As 2 system within 
the energy range -6 eV to 6 eV. PDOS of various orbitals 
like s, p, d and total density of states are indicated by blue, 
green, red solid lines and black dotted line (sum) respectively. 
Vertical dotted line at 0 eV represents the Fermi level. 


II. THEORY AND COMPUTATIONAL 
METHOD 


Our first principles electronic structure calculations 
were accomplished by using plane wave psudopoten- 
tial method based on density functional theory (i^ . 
where the electronic exchange correlation is treated un¬ 
der the generalized gradient approximation (GGA) us¬ 
ing Perdew BurkeErnzerhof (PBE) /PW 91 (d^ func¬ 
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FIG. 2: (Colour online) Calculated normalized total den¬ 
sity of states (red solid line) using VGA approach for 
Bai_a;Ka;Fe 2 As 2 (hole doping) system with (a) x = 0.1, (b) 
X = 0.2, (c) X = 0.3 and (d) x = 0.4 within the energy inter¬ 
val -5.5 eV to 2.5 eV (near Fermi level, indicated by vertical 
dotted line at 0 eV). 



FIG. 3: (Colour online)Calculated normalized total density of 
states (blue solid line) using SC approach for Bai_a:Ka;Fe 2 As 2 
(hole doping) system with (a) x — 0.12, (b) x = 0.19, (c) 
X = 0.31 and (d) x = 0.38 within the energy interval -5.5 eV 
to 2.5 eV (near Fermi level, indicated by vertical dotted line 
at 0 eV). 


tional. A number of authors addressed to the fact that 
density funtional theory within local density approxima¬ 
tion (EDA) as well as generalized gradient approxima¬ 
tion (GGA) was unable to achieve the experimental value 
of ZAs (fractional co-ordinate of As) fairly accurately 
[IE tlE HE HI) EE] • In fact the optimized value of 
ZAs is about 0.1 A smaller than that of the experimental 
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ZAs- The reason behind this discrepancy is the presence 
of strong magnetic fluctuation, emerging due to the exis¬ 
tence of Fe atoms in these materials [5^ . Also it should 
be noted that electronic structures (like Fermi surface), 
calculated using optimized lattice parameters like a, b, 
c and ZAs do not resemble with that of the experimen¬ 
tally measured one [13, [5l|. All these studies force us 
to take experimental lattice parameters i.e., a, b, c and 
ZAs m, Hll, m, [s^l instead of geometry optimized (total 
energy minimization) lattice parameters as the inputs of 
our first principles electronic structure simulations. We 
use room temperature tetragonal experimental lattice pa¬ 
rameters of undoped as well as various doped BaFe 2 As 2 
systems to simulate the electronic structures (band struc¬ 
ture, density of states). Non-spin-polarized single point 
energy calculations were carried out for tetragonal phase 
with space group symmetry I4/mmm (No. 139) using ul- 
trasoft pseudo-potentials and plane-wave basis set with 
energy cut off 500 eV and higher and self-consistent field 
(SCF) tolerance as 10“® eV/atom. Brillouin zone is sam¬ 
pled in the k space within MonkhorstPack scheme and 
grid size for SCF calculation is chosen as per requirement 
of the calculation for various systems and approaches. 
For simulating Fermi surfaces, grid size of SCF calcula¬ 
tion is chosen as 30 x 30 x 41. As mentioned earlier in 
this work we theoretically implement the doping using 
two different approaches (a) Virtual crystal approxima¬ 
tion (VCA) and (b) Super-cell (SC) method. Below we 
briefly discuss the theory, development, advantages and 
disadvantages of the two formalisms. 

Study of disordered systems employing first principles 
electronic structure methods requires some amount of ap¬ 
proximation to implement disorder or doping theoreti¬ 
cally. To treat such disordered or doped systems there 
exist two very popular methods as mentioned earlier: 
super-cell (SC) method and virtual crystal approxima¬ 
tion (VCA). It is very important to mention the advan¬ 
tages and shortcomings of both the approaches before 
we discuss the theoretical background of those meth¬ 
ods. The super-cell method can give more accurate re¬ 
sults than the other one but certainly needs much more 
computational resources as well as time compared to the 
VCA formalism. Accuracy of calculation is associated 
to the fact that the SC method can portray the local 
interaction between two ’real’ atoms whereas the VCA 
method unable to recount that. One of the most impor¬ 
tant drawback for SC method is the artificially imposed 
periodicity which sometime may not mimic the actual 
physical situation of randomly doped system. This ap¬ 
proach certainly . demand for a critical look at the finite 
size effect as well as the imposed artificial periodicity. 
On the other hand, a much lucid and computationally 
inexpensive approach for treating doped systems is to in¬ 
troduce the virtual crystal approximation [53 |. in which 
a doped crystal with the original periodicity, but com¬ 
posed of fictitious virtual atoms is created to mimic the 
actual doped system. To apply VCA method in order to 
study the case of doping, one must concern about certain 


important issues: firstly the accuracy of calculation, sec¬ 
ondly the efficiency for treating the heterovalent systems 
and lastly the electronic property that we are trying to 
calculate. There exists a different formalism called com¬ 
putational alchemy [13, to go beyond VCA. But this 
method is much more intricate than slandered VCA for¬ 
malism, demanding the utilization of density-functional 
linear-response methods. Now we discuss the develop- 



FIG. 4: (Colour online) Calculated band structures (BS) of 
Bai_3;Kj;Fe2As2 (hole doping) system near Fermi level (-2.5 
eV to 2.5 eV) using SC (blue) methods for (a) x = 0.25, (c) 
X = 0.5 and VCA (red) methods for (b) x = 0.3, (d) x = 0.5. 
Horizontal solid black line at 0 eV indicates the Fermi level. 
Same k path (F—A—F—V) has been chosen for BS calculation 
for two different methods. 


ment of SC and VCA formalisms and some of their the¬ 
oretical aspects. The super-cell method defines a new 
unit cell with relatively large number of atoms where 
doping has been implemented, and then (conceptually) 
repeats that cell throughout all space using periodic (or 
toroidal) boundary conditions [56j. For constructing our 
super-cells we follow three steps. At first we replicate 
the crystal unit cell in all 3 directions to get a larger new 
unit cell with relatively large number of atoms. Then we 
introduce doping by substituting atoms with the dopant 
atoms. At the last step, we check the convergence of elec¬ 
tronic properties with gradually increasing the size of the 
super-cell. On the other hand, the base of the VCA for¬ 
malism is simple mixing of pseudopotential The 

supremacy of VCA method is its simplicity but in some 
cases it is not sufficiently accurate. The reason of incor¬ 
rectness lies in the fact that only mixing of pseudopo¬ 
tential is considered. Ramer and Rappe developed a 
more perfect VCA approach by considering averages at 
the atomic level, such as the averages of eigenvalues of 
valence orbitals, charge densities of core electrons , nu¬ 
clear potentials, wavefunctions etc. Major disadvantages 
of this method is that it is incapable of treating heterova¬ 
lent atoms. In this work we use the VCA method, devel- 
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oped by Bellaiche and Vanderbilt based on weighted 
averaging of pseudopotentials. Below we describe the 
method briefly. Total energy of N valence electrons in 
terms of one-particle wave functions (f>i can be written in 
atomic units as: 

Et[{cI),}, {Ri}] = + Ve.t I > 

i 

+ + + ( 1 ) 

where 

(r, r') = ^ Vp',(r - R,, r' - R.) (2) 

I 

R/ is the location of the site 1, and Vp^ are the pseudopo¬ 
tentials. n(r) is the electron density at r, Exc is the 
exchange correlation energy and ?7({Ri}) is the ion-ion 
interaction energy. In case of local as well as non local 
pseudopotential, ’VGA’ operator equation can be derived 
by averaging the pseudopotential of the doped element on 
site 1. 


= (1 +2;Vpf(r,r') (3) 

where, e.g., A=Ba and B=K in Bai_a;K 2 ;Fe 2 As 2 . For 
other undoped sites C (Fe or As site), pseudopotential 
can be written as, 

= ( 4 ) 

Now Vexi takes the form, 

Vextlr, r') = EE - R-is) (5) 

i s 

where and Wg are the pseudopotential and ’weight’ 
(statistical composition) of an atom of type S respec¬ 
tively. 



FIG. 5: (Colour online) Calculated normalized total density 
of states using SC approach (blue solid line) as well as VCA 
approach (red solid line) for (a,b) BaFe 2 (Asi_xPx )2 (isovalent 
doping) system with x = 0.32 within the energy interval -5.5 
eV to 2.5 eV (near Fermi level, indicated by vertical dotted 
line at 0 eV). 



Energy (eV) 


FIG. 6: (Colour online) Calculated normalized total density 
of states for BaFe 2 -a;Coa; As 2 (electron doping) system using 
VCA approach (red solid line) with (a)x = O.f, (b)a; = 0.2 
as well as SC approach (blue solid line) with (a.)x = 0.12, 
(b) 2 ; = 0.25 within the energy interval -5.5 eV to 2.5 eV (near 
Fermi level, indicated by vertical dotted line at 0 eV). 



FIG. 7: (Colour online) Calculated band structures (BS) of 
BaFe 2 -a;Cox As 2 (electron doping) system near Fermi level (- 
2.5 eV to 2.5 eV) using (a) SC (blue)(a; = 0.25) and (b) VCA 
(red) (a; = 0.2) methods. Horizontal solid black line at 0 eV 
indicates the Fermi level. Same k path (F — A — F — V) has 
been chosen for BS calculation for two different methods. 


III. RESULT AND DISCUSSION 

We first calculate the orbital projected density of states 
(DOS) of BaFe 2 As 2 system using experimental tetrago¬ 
nal lattice parameters a, b, c and zas [13 ■ FIG. [1] 

calculated partial density of states (PDOS) as well as to¬ 
tal DOS near Fermi level (-6 eV to 6 eV) are presented. 
Gontributions of various orbitals namely s, p and d to 
the total DOS are depicted by blue, green and red solid 
lines respectively and dotted line indicates the total DOS 
of the system. Fermi energy (E^) is designated by the 
vertical dotted line at 0 eV. It is very clear from FIG. 
[T] that d-orbital of Fe has the major contribution to the 
total DOS near Fermi level. We also observed a signif¬ 
icant contribution of p-orbital near Fermi level, mainly 
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coming from As p-orbital, to the total DOS. These ob¬ 
servations are very much consistent with earlier exper¬ 
imental and theoretical findings [H, As men¬ 

tioned earlier that pure BaFe 2 As 2 system is not super¬ 
conducting at ambient pressure but doping (at any of the 
sites) is one the possible routes to make it a superconduc¬ 
tor. In this work through detail first principles electronic 
structures simulation, we provide a comparative study of 
electronic structures of various doped BaFe 2 As 2 systems, 
where implementation of doping has been done by two 
different theoretical approaches: Virtual crystal approx¬ 
imation (VGA) and Super-cell (SC) methods. We take 



FIG. 8: (Colour online) Calculated normalized total density 
of states for BaFe 2 -a:Ni 2 ,As 2 (electron doping) system using 
VGA approach (red solid line) with {a,)x = 0.1, (b)x = 0.5 
as well as SC approach (blue solid line) with (a)a: = 0.12, 
(b)a; = 0.5 within the energy interval -5.5 eV to 2.5 eV (near 
Fermi level, indicated by vertical dotted line at 0 eV). 

five doped BaFe 2 As 2 systems to examine the evolution of 
DOS and band structure (BS) with doping for both VGA 
and SC methods as mentioned previously. In FIG. [5] and 
FIG. [3] evolution of calculated total normalised DOS of 
Bai_a;Ka;Fe 2 As 2 system near Fermi level as a function 
of doping concentration has been shown for VGA (red) 
and SC (blue) approaches respectively. Calculated band 
structure of Bai_a;Ka;Fe 2 As 2 systems near Fermi energy 
has been shown for x = 0.25 and x = 0.5 using SC 
approach in FIG. 0^ and FIG. 0}: respectively. Same 
has been presented using VGA formalism with a; = 0.3 
and X = 0.5 in FIG. 0 Jd and FIG. 0Ji respectively. BS 
calculation has been performed along same k path for 
both the approaches (VGA and SG) to compare the elec¬ 
tronic structures of these two cases. Calculated BS and 
DOS for Bai_ 2 ,K 3 ;Fe 2 As 2 system reveal that both the 
method (VGA and SC) produce very similar electronic 
structures. Since DOS around Fermi level of BaFe 2 As 2 
system is mainly Fe d-orbital derived, K doping in place 
of Ba does not affect the electronic structure near Fermi 



FIG. 9: (Golour online) Calculated band structures (BS) of 
BaFe 2 -a:NbAs 2 (electron doping) system with x = 0.5 near 
Fermi level (-2.5 eV to 2.5 eV) using (a) SC (blue) and (b) 
VCA (red) methods. Horizontal solid black line at 0 eV indi¬ 
cates the Fermi level. Same k path (F — A — F — V) has been 
chosen for BS calculation for two different methods. 
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FIG. 10: (Golour online) Calculated normalized total den¬ 
sity of states (red solid line) using VCA approach for 
BaFe 2 -a:Rua;As 2 (isovalent doping) system with (a) x = 0.1, 
(b) X = 0.2, (c) X — 0.4, (d) x = 0.6, (e) x = 0.7 and (f) 
X = 0.8 within the energy interval -5.5 eV to 2.5 eV (near 
Fermi level, indicated by vertical dotted line at 0 eV). 


level significantly (with increasing K doping concentra¬ 
tion FSs become more two dimensional [^). As a result, 
K doping only shifts the chemical potential indicating the 
appearance of extra hole to the system which is reflected 
in the calculated as well as experimentally measured FSs 
But DOS profile as well as BS is almost 
unaffected by K doping, in contrast to the cases of dop¬ 
ing at other sites. This is the main reason of getting 
same electronic structures for VCA and SC approaches. 
Central message from this study is that instead of SC 
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method one can use time saving and numerically simpler 
VGA method for studying the electronic structures such 
as BS, DOS, Fermi surface etc. of hole doped BaFe 2 As 2 
systems. Next we have studied the case of isovalent P 
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FIG. 11: (Colour online) Calculated normalized total den¬ 
sity of states (blue solid line) using SC approach for 
BaFe 2 -a:Rua;As 2 (isovalent doping) system with (a) a: = 0, 
(b) X — 0.25, (c) X = 0.5, (d) x = 0.75, (e) a; = 1 and 
(f) a: = 2 within the energy interval -5.5 eV to 2.5 eV (near 
Fermi level, indicated by vertical dotted line at 0 eV). 

doping in place of As for BaFe 2 As 2 system. In FIG. [5^ 
and FIG. [5b we present the normalized DOS of 32% P 
doped BaFe 2 As 2 systems for SG (blue) as well as VGA 
(red) approaches respectively. Apart from nominal con¬ 
trast in DOS profile very near to the Fermi level, we 
observed hardly any difference between the DOS calcu¬ 
lated from SG and VGA methods. As with P doping 
in As site effectively modifies the Fe-As hybridization, 
one should expect some moderation of DOS around E^r, 
but unaltered large Fe-d orbital contribution to the total 
DOS restricts any notable modification of total DOS near 
Fermi level. As far as the two distinct methodologies are 
concerned, we found both the approaches produce very 
similar DOS profile within the energy interval -5.5 eV to 
2.5 eV for 32% P doped BaFe 2 As 2 systems. As discussed 
earlier DOS near EiT- is dominated by Fe d-orbital, thus 
methods of implementation (VGA and SG) of P doping 
in As site play insignificant role in the calculation of total 
DOS. Later we will discuss this issue to shed more light 
on the differences in the calculated electronic structures 
by these two approaches, providing the results of calcu¬ 
lated PDOS of the 32% P doped BaFe 2 As 2 system. Now 
we study the case of doping on Fe site which is very inter¬ 


FIG. 12: (Colour online) Calculated band structures (BS) of 
BaFe 2 -a:Rua;As 2 (isovalent doping) system near Fermi level (- 
2.5 eV to 2.5 eV) using (a,c) SC (blue) method with x = 0.25 
and X = 0.5, (b,d) VGA method (red) with x = 0.2 and 
X = 0.5. Horizontal solid black line at 0 eV indicates the 
Fermi level. Same k path (F — A — F — V) has been chosen 
for BS calculation for two different methods. 



FIG. 13: (Colour online) Calculated d (red) and p (green) 
projected partial density of states (PDOS) near Fermi level 
(-5.5 eV to 2.5 eV) for BaFe 2 (Asi-a;Pa :)2 (isovalent doping) 
system with x = 0.32 using (a) SC and (b) VGA approaches. 


esting as electronic structure near the Fermi level is dom¬ 
inated by Fe-d orbitals. First, we study the case of Go 
and Ni doping in place of Fe which exerts additional elec¬ 
tron to the system and chemical potential shifts accord¬ 
ingly, is referred as electron doping. Phase diagram of 
electron doped (Go/Ni) Fe-based superconductors clearly 
suggests the presence of superconducting as well as var¬ 
ious other phases like SDW, Nematic etc. It should be 
noted that in case of electron doping, superconducting 
regime (also other phases) extended upto 15 % Go/Ni 
doping concentration. So study of electronic structure 
in the low doping regime is well enough to draw conclu¬ 
sion about superconducting as well as other exotic normal 
state properties involving electronic structures. For elec¬ 
tron doped system, we have presented two case studies: 
(1) BaFe 2 -a:Goa:As 2 systems and (2) BaFe 2 -a:Nia;As 2 sys- 
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terns. In FIG. [6^ and FIG. IUd, we display the simulated 
normalized DOS of BaFe 2 _a;Go 3 ;As 2 systems with x = O.I 
and X = 0.2 and in FIG. [8^ and FIG. Eb, we depict the 
calculated normalized DOS of BaFe 2 -a;Nia;As 2 systems 
with X = Q.l and x = 0.5, where doping has been im¬ 
plemented through VGA method. FIG. Et and FIG. [6ji 
manifest the normalized DOS for BaFe 2 -a;CoxAs 2 sys¬ 
tems with X = 0.12 and x = 0.25, calculated using SC 
approach. Same has been presented for BaFe 2 -xNia;As 2 
systems with x = 0.12 and x = 0.5 in FIG. Et and 
FIG. EJl- In FIG. [7] we demonstrate the calculated BS 
of BaFe 2 -a;Coa;As 2 systems with x = 0.2 for VGA (red) 
and with x = 0.25 for SC (blue) methods. FIG. [9b and 
FIG. Eb represents the simulated BS of BaFe 2 -xNia;As 2 
systems with x = 0.5 for SC and VGA approaches re¬ 
spectively. It is evident from calculated BS as well as 
calculated DOS (FIG. |6b, FIG. Hi, FIG. |Zl FIG. H) that 
VGA and SC approaches produce very similar electronic 
structure near Ep’ in the low doping regime but dissimi¬ 
larities in the electronic structure arises in case of large 
doping like 50 % Ni doped BaFe 2 As 2 system. In this 
case calculated DOS using VGA method has four distinct 
peaks within the energy interval -5.5 eV to 2.5 eV (see 
fig FIG.Eb)- Where as calculated DOS implementing SC 
approach has five peaks (see fig FIG.Eb) and DOS profile 
very near to the Fermi level is completely different from 
that of the calculated by VGA method. In SC approach 
50 % doping has been accomplished by substituting 50 
% of the Fe atoms by Ni atoms from various random lat¬ 
tice sites of the bigger unit cell that we call super-cell. 
Consideration of two distinct sites for Ni and Fe in SC 
method is the main reason of the occurrence of two peak 
near Fermi energy instead of one. In VGA method 50 % 
doping has been implemented by mixing Fe and Ni pseu¬ 
dopotential in equal weightage. That is why one gets one 
peaks in the DOS profile in VGA approach at a partic¬ 
ular energy as a result of averaged peudopotential of Ni 
and Fe. Later we will show through PDOS calculation, 
in SC method Ni and Fe peaks appear at slightly differ¬ 
ent energies on contrary to the case of calculated DOS by 
VGA method. Calculated BS (50 % Ni doped BaFe 2 As 2 
system) for the two approaches appeared to be very sim¬ 
ilar as far as position of bands are concern but a large no 
of extra band arises due to implementation of Ni atom in 
distinct lattice sites in SC calculation (see FIG. H). 

Next we study the case of isovalent doping on Fe site 
i.e., Ru substitution in place of Fe. FIG.ITOldemonstrates 
the calculated normalised DOS of Ru doped BaFe 2 As 2 
system as a function of Ru concentration where doping 
has been implemented via VGA formalism. On the other 
hand, in FIG. |TT]we have presented the simulated nor¬ 
malized DOS for various Ru concentrations considering 
SC approach for the execution of Ru doping. FIG. ITT] de- 
picts the band structures of Ru doped BaFe 2 As 2 system 
for two different Ru doping concentrations using both the 
two formalisms: VGA (red) and SC (blue). From FIG.fTTl 
one can see that the position of the bands are not same 
for two different formalisms. It is very clear from these 
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FIG. 14: (Colour online) Calculated partial density of 
states (PDOS) near Fermi level (-5.5 eV to 2.5 eV) for 
BaFe 2 -a;Cox As 2 (electron doping) system using SC method 
with (a) x = 0.12 and (b) x = 0.25 indicating contribution of 
Co-d (red) and Fe-d (orange) orbitals as well as total d-orbital 
(green) to total DOS and VC A method with (c) x = 0.1 and 
(d) X = 0.2 indicating PDOS of d-orbital (blue). 


figures (FIG. [TOl FIG. [TTl FIG. [12]) that the simulated 
DOS and BS are very much different for VGA and SC 
formalisms. It appears that in case of VGA approach the 
chemical potential shifts a bit for larger doping concen¬ 
tration and deviates from the electronic structures cal¬ 
culated using SC method. A more detailed inspection of 
electronic structure is thus required to unfold the reason 
of dissimilarities between the electronic structures calcu¬ 
lated by VGA and SC formalism. Below we present our 
results of calculated PDOS of all the BaFe 2 As 2 systems 
that we discussed earlier to get a better insight of the 
electronic structures calculated by these two approaches. 

Since K doping at Ba site does not alter the DOS 
much, near the Fermi level, it is irrelevant to show the 
PDOS of these K doped systems or more precisely hole 
doped systems in the context of our main aim of this 
work. First we see the differences in calculated PDOS 
for VGA as well as SC methods in case of 32% P doped 
BaFe 2 As 2 system. FIG. [13] depicts the calculated orbital 
projected PDOS for the two approaches. Total d and p 
orbital contributions are presented for both the cases. 
P doping in As site effectively modifies the PDOS of p- 
orbital which in turn results some amount of moderation 
of Fe-As hybridization. From FIG. [13] it is evident that 
both the approaches simulate same PDOS of d-orbital 
but there is remarkable differences in PDOS of p-orbital 
calculated by two different methods. The differences 
are more pronounced if one compares the PDOS of 
p-orbital with respect to the PDOS of d-orbital for the 
two formalism, that we considered. Now we proceed to 
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FIG. 15: (Colour online) Calculated partial density of 
states (PDOS) near Fermi level (-5.5 eV to 2.5 eV) for 
BaFe 2 -a;Nia;As 2 (electron doping) system using SC method 
with (a) X = 0.12 and (b)x = 0.5indicating contribution of 
Ni-d (red) and Fe-d (orange) orbitals as well as total d-orbital 
(green) to total DOS and VCA method (c) x = 0.1 and (d) 
X = 0.5 indicating PDOS of d-orbital (blue). 


the case of Co/Ni doping in Fe site. In FIG. [Ml and FIG. 
[T5] we depict the calculated PDOS of BaFe 2 -a;Goa^As 2 
and BaFe 2 -xNixAs 2 systems respectively. In the left 
column of both the Figures, results of SC calculation is 
presented, where PDOS of Fe-d (orange), Ni/Co-d (red) 
as well as total d-orbitals (solid green line) are shown 
for various doping concentrations as indicated in the 
figures. Calculated total d-orbital contribution (solid 
blue line) using VCA method is also displayed in the 
right column of both the figures (see FIG. [14] and FIG. 
dS]) for comparison with that of the calculated PDOS 
using SC methods (left column). Comparative study of 
calculated PDOS of Co/Ni doped systems for both the 
methods reveals that in case of 50% Ni doping VCA 
and SC method produce significantly different electronic 
structures. One can see from FIG. dSiD that the peak 
of the d-orbitals of Ni and Fe are at slightly different 
energies and sum of these two DOS profiles produce the 
total d-orbital contribution. But in case of VCA, PDOS 
of the total d-orbital as shown in FIG. dSH is totally 
different from that of the SC calculation. Although 
implementation of higher doping is not that much 
important in case of electron doped system as discussed 
earlier but one must concern about the fact that in case 
of higher doping concentration VCA method is unable 
to simulate the accurate electronic structures. Next 
we show in FIG. the simulated PDOS of Ru doped 
BaFe 2 As 2 system for VCA and SC methods. In SC 
calculation (left) we presented the d-orbital contribution 
of Ru-d (red) as well as Fe-d (orange) separately and 
total d orbital contribution displayed by solid green 
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FIG. 16: (Golour online) (a,b,c) Calculated partial density 
of states for BaFe 2 -xRuxAs 2 (isovalent doping) system with 
X = 0.25(a), X = 0.5(b) and x = 0.75(c) using SC method. 
Contribution of Fe-d (orange), Ru-d (red) and total d (green) 
orbitals are indicated. (d,e,f) Normalized PDOS of d-orbital 
(blue) for BaFe 2 -a:Rua;As 2 (isovalent doping) system with x = 
0.2(d), x = 0.5(e) and x = 0.8(f) using VCA method. 


line (see FIG. [TBh . b, c). On the right side of FIG. [ini 
total d-orbital contribution has been depicted for VCA 
method. Comparison of both the electronic structure 
reveals that in VCA formalism, mixed d-orbital of Ru 
and Fe atom deviates from the PDOS of total d-orbital 
calculated by SC method. Although Ru doping at the 
Fe site is iso-electronic substitution but total number of 
electron as well as size of Ru is much greater than that 
of Fe. Furthermore Ru has 4d orbital as the outermost 
orbital, whereas outermost orbital of Fe is 3d. It also 
appears that in the low doping concentration, apart 
from an additional shift of the chemical potential in 
the calculated electronic structures by VCA method, 
there is not much difference in the electronic structures 
calculated by these two formalisms. But in case of higher 
doping concentrations, electronic structure simulated 
using VCA formalism deviates remarkably from that 
of the SC method. Perhaps in VCA formalism, the 
mixing of 3d (Fe) orbitals with 4d (Ru) orbitals which 
are more broader than 3d, is not implemented properly 
resulting an additional shift in the chemical potential 
in DOS profile for VCA approach. This discrepancy 
can be overcome in the low doping regime. In FIG. 
[TZl (left column) we depict the calculated FSs of Ru 
doped BaFe 2 As 2 systems using VCA formalism. In 
the right column, we present the FSs of the same by 
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FIG. 17: (Colour online) Simulated Fermi surfaces of 
BaFe 2 -a:Rua;As 2 using VGA methods with out any chemical 
shift (left column) and with proper chemical shift (right col¬ 
umn). Chemical shift is estimated by comparing calculated 
DOS for VGA and SC approaches. 


giving an additional shift , which is estimated from 
the comparison of calculated DOS by VGA and SC 
methods. It is clear from FIG. [m that in the low 
doping regime the shifting of chemical potential work 
very nicely to produce the FSs, that resembles with the 
experimentally measured FSs as well as calculated FSs 
using SC methods [5l|, |H, . But in case of higher Ru 

doping concentration this shifting of chemical potential 
does not work to reproduce the FSs accurately as 
evident from the simulated FS in the last row of FIG. nzi 


Calculated FSs using VGA formalism in case of 40% Ru 
doping (last row of FIG. ITTl) certainly not resembles with 
that of the SC calculation on 50% Ru doped systems [5l|. 

IV. CONCLUSION 


Implementing VGA and SC methods for doping, we 
have presented a detail comparative study of electronic 
structures of various doped BaFe 2 As 2 systems through 
first principles simulation. We have investigated the den¬ 
sity of states (including partial density of states), band 
structures and Fermi surfaces of electron doped, hole 
doped as well as isovalently doped BaFe 2 As 2 systems 
for VGA and SC approaches. In case of hole doping 
and iso-electronic P substitution on As site, both the 
methods produce very similar electronic structures as far 
as DOS and BS calculations are concerned. Calculated 
PDOS of this system give clear indication about the fact 
that P doping in As site (isovalent) moderate Fe-As hy¬ 
bridization. So in this case there is some difference in 
the calculated PDOS for VGA and SC approaches. In 
case of electron doped (Co/Ni) BaFe 2 As 2 systems VGA 
and SC methods provide very similar electronic struc¬ 
tures for both the cases for low doping concentrations. 
But in case of higher doping concentration, we find that 
the calculated electronic structures are remarkably differ¬ 
ent for these two methods as clearly evident from the BS 
and PDOS calculation of 50% Ni doped BaFe 2 As 2 sys¬ 
tems. On the other hand, in case of iso-electronic substi¬ 
tution in Fe site (Ru substitution in Fe site), calculated 
electronic structure is very much dissimilar for different 
methods. An extra shift in the chemical potential is ob¬ 
served in the electronic structures of Ru doped BaFe 2 As 2 
systems, when VGA method is employed for doping. But 
for higher Ru doping concentration the electronic struc¬ 
tures are quite different for these two methods. We also 
explicitly show that by considering this chemical shift in 
a proper way, one can simulate FSs accurately by VGA 
formalism in case of low Ru doping concentration. This 
work give a clear idea about the application of VGA and 
SG methods for calculating the electronic structure of 
various doped 122 systems of Fe-based superconductors. 
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